Theoretical study of magnetic properties and hyperfine interactions in cr-FeV alloys 
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Electronic structure Korringa-Kohn-Rostoker calculations for the cr-phase in Feioo-zV^ were per- 
formed in compositional range of its occurance (~ 34 < x <~ 65). Fe and V magnetic moments 
and hyperfine fields were determined for five inequivalent lattice sites in two models of magnetic 
structure, namely ferromagnetic FM and so-called anti-parallel one, APM, dominated by antifer- 
romagnetic coupling. The average magnetic moments calculated for FM state overestimate the 
experimental data, whereas the corresponding quantities computed for APM state, underestimate 
them. Such a behavior remains in line with total energy values being similar for both models. The 
calculations showed that both average magnetic moments and hyperfine fields (on Fe and V atoms) 
vary with a number of Fe atoms in the nearest neighbor shell, NNfb, starting from critical values 
of NNF e , characteristic of each site. 

The calculated hyperfine fields for Fe and V showed an important role of valence contributions, 
being strongly dependent on local magnetic moments arrangements. In the case of Fe, the computed 
hyperfine fields are overestimated with respect to the Mossbauer data, while the corresponding values 
for V are underestimated with respect to the NMR results. However, the linear correlation between 
average magnetic moment and hyperfine fields, observed experimentally in FeV cr-phase, can be 
well reproduced when combining theoretical results for the two above-mentioned magnetic structure 
models. 



I. INTRODUCTION 



Among the so-called Frank-Kasper phases that can be 
geometrically described in terms of the basic coordina- 
tion polyhedron with coordination number equal to 12, 
14, 15 and 16^, the cr-phase that can occur in transition- 
metal alloy systems, is known as the one without definite 
stoichiometric composition. Consequently, it can exist in 
different composition ranges and, therefore, its physical 
properties can be tailored by changing its constitutional 
elements, and within a given constitution, by changing 
chemical composition. Among about 50 examples of bi- 
nary cr-phases that have been reported so far in the liter- 
ature, only two viz. er-FeCr and cr-FeV are known to have 
well evidenced magnetic properties^—, though a knowl- 
edge of the cr-phase magnetism remains not complete. In 
particular, neither values of the magnetic moments, fi, 
on particular lattice sites nor a coupling between them is 
known. A lack of stoichiometry, leading to a huge number 
of different atomic configurations, combined with a weak 
magnetism i.e. small values of fj, makes it very difficult 
for microscopic methods like Mossbauer spectroscopy or 
nuclear magnetic resonance to uniquely determine spin- 
densities at the level of the sublattices. Non-availability 
of big enough single-crystals of the cr-phase prevents the 
neutron diffraction techniques to be used for that pur- 
pose. In these circumstances, performing theoretical cal- 
culations aimed at determining electronic and magnetic 
structures of the phase is highly needed and helpful. 

The present paper reports such calculations for the 
Feioo-icVa; system where the cr-phase exists in a wide 
range of composition (~ 34 < x <~ 65), and its mag- 
netic properties show a strong x- dependence. For exam- 
ple, the Curie temperature may vary between K, for 



x « 65, and ^320 K, for x w 34 8 . This gives a unique 
chance for testing the calculations which were already 
successfully applied for similar purpose in the cr-FeCr al- 
loy system^ where the phase exists in a much narrower 
compositional and temperature ranges. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations for cr— FeV com- 
pounds have been performed using the charge and spin 
self-consistent Korringa-Kohn-Rostoker (KKR) method. 
The muffin-tin crystal potential has been employed 
within the local density approximation (LDA) frame- 
work, where the von Barth-Hedin formula^ for the 
exchange-correlation part was used. The self-consistent 
crystal potentials were converged below 0.1 mRy and 
charges below 10 -3 e. The spin-polarized densities of 
states were computed employing the tetrahedron k-space 
integration technique. The local magnetic moments (of 
Fe and V) inside each muffin-tin spheres were determined. 
Accordingly, the Fermi contact hyperfine fields (Hhf) re- 
sulted from the spin densities extrapolated to r — 0, p(0), 
with the use of the well-known non-relativistic formula 
H hf = 3^5/9(0) « 524p(0) in kOeU. The KKR cal- 
culations allowed to determine core (H core ) and valence 
(H va i — Hhf — H core ) terms, neglecting however dipo- 
lar and orbital contributions to the total hyperfine field. 
It is worth noting that electronic structure calculations 
have been done for ordered approximants of the disor- 
dered cr— FeV alloys, and for that purpose the symmetry 
of the unit cell was lowered to a simple tetragonal one. In 
practice, the tetragonal unit cell (space group P42/nmm) 
and atomic positions were unchanged but variable occu- 
pancy made all 30 atomic positions crystallographically 
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FIG. 1: (Online color) Magnetic moments of Fe and Cr atoms for five crystallographic sites versus the number of NNFe atoms 
as obtained with the FM model of the magnetic ordering. Solid lines stand for the best fits to the data. 
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FIG. 2: (Online color) Magnetic moments of Fe and Cr atoms for five crystallographic sites versus the number of NNFe atoms 
as obtained with the APM model of the magnetic ordering. Solid lines stand for the best fits to the data. 



nonequivalent (space group PI) being next occupied ex- 
clusively either by Fe or V atom. Calculations have been 
performed for two er-FeV alloys with different composi- 
tions, namely Fe2oVio and Fe^Vis. They are equivalent 
to the formula Feioo-xV x with x = 33.3 and x = 60.0, 
respectively. The chosen compositions correspond fairly 
well to border concentrations of the cr-phase existance in 
the Fe-V alloy system^. Fe and V atoms were distributed 
over five inequivalent sublattices of the tetragonal unit 
cell in close agreement with experimentally determined 
occupancies^. Actually, the used here atomic configura- 
tions in the unit cell remained strictly the same as those 



considered previously in electronic structure calculations 
of the cr-FeV phase in a paramagnetic stated. 



III. RESULTS AND DISCUSSION 

A. Magnetic Moments 

In order to study magnetic properties from first prin- 
ciples calculations, two collinear models of a possible 
magnetic ordering were taken into account: (i) parallel 
alignment of local magnetic moments, called here fer- 



3 



romagnetic, FM, as derived directly from spin-polarized 
calculations, and (ii) antiparallel alignment, called here 
APM. In the latter, Fe/V atoms in selected sublattices 
were initially coupled antiferromagnetically according to 
predictions of a symmetry analysis as given in details 
elsewhere 9 . The magnetic moments and hyperfmc fields, 
Hh f , corresponding to both Fe and V atoms distributed 
over the five sublattices were a subject of a further anal- 
ysis. Their non-zero values were found for both Fe and V 
atoms. As a rule, the magnetic moments (and the hyper- 
fine fields) on V atoms were found to be much smaller, 
and, generally, polarized opposite, to those on Fe ones. 

The calculated values of fi and \Hhf\ for each sublat- 
tice were analyzed versus a number of Fe atoms in the 
nearest neighborhood, NN Fe . Such approach combining 
the results obtained for the two border concentrations, 
shows that the two sets of the computed data are com- 
plementary, and they cover the expected range of NN Fe 
for all sublattices. 

A detailed analysis of the fi(N Np e )-dependence of the 
magnetic quantities for each of the sublattices shows that 
this relationship can be split, in most cases, into two 
regions showing a linear character. In the investigated tr- 
FeV system, an average magnetic moment per an atom, 
< fi > , on a given sublattice remains zero up to a critical 
value of NN Fe = NN crit , and for NN Fe > NN crit it 
varies linearly in accord with the equation 

f a(NN Fe - NN crlt ) for NN Fe > NN crlt 
M \ for NN Fe < NN crlt 

(1) 

A similiar dependence on NN Fe was also obtained for 
the average \Hhf\, < \Hhf\ > (NN Fe ), for all sublat- 
tices, and both considered models of the magnetic order- 
ing (FM and APM). 

Noteworthy, a similar discontinous [i(NN Fe )- 
dependence was already reported in the analysis of the 
atomic moment in bcc Fe-V alloys^, and also used in 
calculations of the hyperfine field distribution curves 
for nano crystalline bcc Fe-V alloys^. The critical value 
of NN Fe as found by these authors was rather small, 
NN cr it = 3, when accounting for the magnetic moments 
localized on Fe-atoms. In the present case of the er-phase, 
NN crit values and slopes a of the lines were found by a 
fitting the data calculated for Fe and V atoms located 
on five various sublattices, separatelly. The obtained 
N N cr it-vahies in that way do not exceed 6, as can be 
clearly seen in Figs. [T] and [5] It is worth noting that 
negative values of NN cr it, as found in several cases, 
mean that the magnetic moments do not vanish for all 
A^Fe-values, as derived from Eq. 1. 

The calculated magnetic moments for both types of 
atoms located on all crystallographic positions combined 
with calculated probabilities of these position occupan- 
cies allowed to determine the average values of the mag- 
netic moments, < fi >, for each sublattice and for any 
concentration of vanadium. In the same way, one could 
compute the average magnetic moment per atom (per 
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FIG. 3: (Online color) Average magnetic moments as calcu- 
lated in FM (triangles) and APM (circles) models, as well as 
experimentally determined data (squares)*. Solid lines are for 
guide eye only. 



unit cell). In practice, such calculations were performed 
for four selected concentrations that corresponded to the 
samples on which the magnetic measurements were per- 
formed. A comparison between the calculated and exper- 
imentally determined quantities is displayed in Fig. [3] 
It can be clearly seen that the values of < fi > calcu- 
lated applying the FM model are systematically larger 
than those found experimentally, < [i exp >. The dif- 
ference between the former and the latter changes from 
about ~ 0.35yUs per atom for x = 34.4 to ~ 0.25/xs for 
x = 55.0. However, the overall dependence of < fi > 
on the vanadium content follows well that found experi- 
mentally. On the other hand, the < /i >-values derived 
from the APM model are underestimated in comparison 
to those of < iJ. eX p > but the agreement bewteen theo- 
retical and experimental values markedly improves with 
increasing V content. For the vanadium richest sample, 
Fe45V55, the measured value of < [i eX p >= 0.05/is is 
practically the same as the one calculated with the APM 
model, < fi >= O.Oi^B- 

The relationship between the measured and calculated 
values of < \i > can be explained in terms of a total 
energy of the cr-phase, E, calculated for the two models 
of magnetic of ordering (E F m and Eapm)- The energy 
difference AE — Eapm — E F m, was determined using 
the electronic structure calculations for the two border 
concentrations of the cr-phase within both FM and APM 
models. In the case of the c-Fe2oVio alloy Ai5-value 
ranges from ~ 0.7 to ~ 1.5 mRy per atom, that would 
favour the FM model. For the cr-Fe^Vig alloy the AE- 
value almost vanishes (less than 0.05 mRy). The total en- 
ergy difference of the cr-phase can be compared to KKR 
results obtained for FM and (hypothetical) AFM order- 
ing in a simple a-phase with the same Fe-V composition. 
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FIG. 4: (Online color) \Hy\ and |i?p e |-values for five crystallographic sites versus the number of NNFe atoms as obtained with 
the FM model of the magnetic ordering. Solid lines stand for the best fits to the data. 
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FIG. 5: (Online color) \Hv\ and |//Fe|-values for five crystallographic sites versus the number of NNFe atoms as obtained with 
the APM model of the magnetic ordering. Solid lines stand for the best fits to the data. 



It yields the value of ~ 15 mRy per atom (also favoring 
FM ordering), being thus an order of magnitude larger 
than in the former case. 

Hence, the small difference between total energies com- 
puted for the FM and APM structures may suggest that 
some regions of the sample have the antiparallel magnetic 
order (APM), whereas in others the magnetic moments 
arrange mostly ferromagnetically (FM). In other words, 
the magnetic properties of the cr-FeV can be thought of 
as a coexistence of two types of the ordering. The rel- 
ative contribution of them should depend on the alloy 
composition. In view of the total energy analysis, the 



contribution of the APM-like regions should be increas- 
ing with the vanadium concentration, x, since the AE 
value decreases with x. The latter effect is also observed 
experimentally in the dependence < ^ exp > (x) - Fig. [3J 

On the other hand, it should be remembered that both 
FM and APM models do not strictly refer to the ferro- 
magnetic and antiferromagnetic ordering of Fe/V mag- 
netic moments in the cr-phase. They should be rather 
treated as starting configurations, and finally converged 
magnetic moment values and their orderings may differ 
from the initial alignments. In consequence, magnetic 
moments on Fe/V atoms may vary both in magnitude 
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FIG. 6: (Online color) (a) Average absolute values of hyper- 
fine fields on Fe atoms as calculated in FM (triangles) and 
APM (circles) models, as well as experimentally determined 
data (squares)^. Solid lines are to guide the eye only, (b) Av- 
erage absolute values of hyperfine fields on vanadium atoms 
for five crystallographic sites as calculated in FM (dashed 
lines) and APM (solid lines) models, as well as experimen- 
tally determined data (circles)^ 7 -. 



and in a mutual coupling, also due to the fact that the 
cr-FeV phase is a chemically highly disordered system. In 
particular, it is possible that the real magnetic ordering 
is not collinear. A necessary condition for the occurence 
of spin-canting i.e. competition of ferromagnetic and an- 
tiferromagnetic interactions, has been here, and also in 
the case of the a—FeCv^ theoretically suggested to exist. 



B. Fe and V hyperfine fields 

Another quantity pertinent to magnetic materials that 
has been calculated is a hypcrhnc held, Hhf. Although 
its correlation with the magnetic moment is usually non- 
linear and, in some instances, even non-monotonous^, 
in the case of the cr-FeV alloys \Hhf\ and fj, exhibit a 
perfect linear correlation^. In addition, for the cr-FeV 



compounds the hyperfine fields had been measured both 
on 57 Fe and 51 V— so they can be used as a proper data 
set for testing pertinent theoretical calculations. Those 
obtained presently by means of the KKR calculations for 
the FM and APM models are shown in Fig. [4] and [5j 
respectively. As can be clearly seen one can distinguish 
two rangies for all sublattices and occupations, namely 
zero- field range (Hhf = 0) for small values of NNp e , and 
proportionality range (H ~ NNp e ) for higher NNp e - 
values. In order to compare the calculated rT^-values 
with the experimental ones (Mossbauer and NMR data) 
we consider only the absolute values, \Hhf\- 

Concerning the Fe-site hyperfine fields, we can make 
a comparison - see Fig. |6K - only between their aver- 
age values, < \Hp e \ >, due to the fact that experimen- 
tally only such quantity is available from the 57 Fe-site 
Mossbauer spectroscopic measurements 8 . As it follows 
from Fig. [6K, the average hyperfine field on Fe atoms as 
calculated in the frame of the FM model decreases lin- 
early with x from < \Hp e \ >~ 155 kOe for x = 34.4 
to < \Hp e \ >~ 60 kOe for x — 55. The corresponding 
< \Hp e \ >-values obtained by means of the APM model 
also decrease linearly with x from < \Hp e \ >« 105 kOe 
to < \Hp e \ >~ 40 kOe, but at a slightly lower rate. The 
measured values of < |izFe| > {%) show a non-linear be- 
havior. For the lowest V concentrated sample (x — 34.4) 
its value lies between those found for the FM and APM 
models. For x — 40 it is equal to the value calculated 
with the APM approach, while for higher ^-values it is 
smaller than the calculated values. The observed quan- 
titative agreement between the measured and the calcu- 
lated values of the hyperfine fields is not fully satisfying 
what can have different reasons. First of all, we should 
bear in mind that in our present calculation only the 
Fermi contact term was taken into account, while other 
possible contributions like dipolar and orbital ones, were 
neglected. But we do expect that the most important 
influence on the computed ii), /-values could have a lo- 
cal ordering of magnetic moments which mostly modify 
its valence terms. It is well seen on the dependencies of 
Hcore and H va i contributions vs. magnetic moments for 
FM model (Fig. [7]). As expected, H core is perfectly linear 
vs. /i with the expected slope of about 10 T//j,b (both 
for Fe and V atoms whatever the sublattice), while for 
H va i there is no apparent correlation with the magnetic 
moment, since it strongly depends on atomic and mag- 
netic surroundings. In a real magnetic structure, valence 
contributions to f may vary both in magnitude and in 
direction. The latter is not provided by the present cal- 
culations when using simplified FM and APM models. 
It is worth noting that for Fe atoms the calculated Hhf 
are dominated by H core (best seen on A, D and E sites 
in Fig.[7|), which are fully correlated with local magnetic 
moments. Since the calculated hyperfine fields overes- 
timate the experimental data, one should expect that 
the valence contributions, which are opposite in sign to 
the core contributions, may diminish this disagreement. 
Hence, accounting for a more complex magnetic struc- 
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FIG. 7: (Online color) H CO r& and i-/„ a i-values for five crystallographic sites versus the magnetic moment fi as obtained 

with the FM model in cr-FeV phases. 
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FIG. 8: (Online color) Average hyperfine field on Fe-atoms, 
< \Hfb\ >, versus average magnetic moment per atom, fi, 
in the cr-FeV phase, as measured and calculated in FM and 
APM models. 



ture (with different local magnetic moments directions 
that mostly affect H va i) could improve experiment ver- 
sus theory agreement for the results on Fe atoms. 

On the other hand the correlation between the calcu- 
lated average Fe-site hyperfine field and the correspond- 
ing magnetic moment (per atom) - see Fig. [8] - agrees 
quite well with the measured one: the slope of the ex- 
perimentally found < \Hp e \ > — < n eX p > line has the 
average value of ~ 19 T//zb if the magnetic moment is 



calculated per an atom in the unit cell (or 14.3 T/fJ-B 
when the magnetic moment is determined per one Fe- 
atom in the unit cell^), while that of the theoretically 
determined lines amounts to ~14 T/hb for the FM or- 
dering and to ^27 T/^b for the APM one. In other 
words, the experimentally found < \Hp e \ > — < ji exp > 
line lies between the two theoretically calculated ones. 
This can be taken as a further support for our approach 
that the magnetic ordering in the studied system can be 
treated as a combination of FM- and APM-ordered re- 
gions. 



Concerning the 51 V fields, we have calculated the fields 
for four sublattices (B, C, D, and E) as we can compare 
these quantities with experimentally determined ones for 
each of the five sublattices^. In addition, the average 
values of the V-site field, < \Hy \ >, were evaluated, too. 
A comparison between the theoretical and experimental 
data, displayed in Fig. gives an evidence that the 
fields calculated for particular sites are underestimated 
by a factor of ~ 2. In consequence, the average field, 
namely the calculated < \Hy\ > decreases with x at 
a lower rate than the measured one. Using the same 
arguments as for the Fe-site hyperfine fields, we expect 
that H va i are not sufficiently well evaluated within the 
FM and the APM models. Fig. [7] shows that absolute 
values of H va i are slightly larger than the corresponding 
ones of H core , resulting in relatively small total hyperfine 
fields on V. Consequently, the calculated hyperfine fields 
became underetimated with respect to the experimental 
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FIG. 9: (Online color) 57 Fe Mossbauer spectra recorded on 
the (a) cr-Fe65.6V34.4 and (b) <r-Fe6o.iV39.9 sample at 4.2 K. 
Dashed lines represent best fits to the data calculated within 
the FM and the APM models. 

C. Mossbauer Spectra 

The calculations reported in this paper can be also 
used in the analysis of Mossbauer spectra. Towards this 
end we assumed the Mossbauer spectrum to be composed 
of five subspectra correspondig to the 57 Fe atoms located 
on the five inequivalent lattice sites. The relative area of 
each subspectrum should be proportional to the num- 
ber of Fe atoms occupying a given site and it was as- 
sumed to be equal to that known from neutron diffrac- 
tion experiments^. Each subspectrum represents a dis- 
tribution of the hyperfme fields (hfd) with quadrupole 
splittings, QS, and relative isomer shift- values, IS, iden- 
tical to those determined from room temperature mea- 
surements and presented in Ref. [Til The hyperfme 
field distributions for each concentration were computed 
on the basis of probability distributions and calculated 
< \Hp e \ > (NNp e ) dependences obtained from FM 
and APM models. As can be clearly seen in Fig. |H1 
the Mossbauer spectrum recorded on the sample with 
x = 34.4 could be cuccessfuly fitted assuming the hfd to 



be weighted average of hfd obtained from the FM model 
(weigt 58%) and the APM one (42%). For the sample 
with x — 40 it was enough to take into accout the hfd 
calculated for the APM model, only. For more concen- 
trated samples it was not possible to get any reasonably 
good fit on the basis of the calculated values. 



IV. CONCLUSIONS 

The electronic structure calculations of ordered ap- 
proximants of the tr-FeV phase were done for border 
compositions of this phase, namely Fe2oVio and Fe^Vis- 
The magnetic properties, i.e. total and atomic magnetic 
moments, as well as Fe and V hyperfme fields were calcu- 
lated using two models for the magnetic moments align- 
ment and compared to experimental values. It was found 
that both fj,(NNpe) and \Hhf \ are zero below NN cr u, 
characteristic of each lattice site, and increase more or 
less linearly for higher NNp e - values. The calculated 
< fi >-values overestimate the experimental data when 
using the FM model, and underestimate them in the case 
of the APM one. This discrepancy was discussed in terms 
of the total energy difference as calculated for the FM 
and the APM models. The analysis of the calculated 
Fermi contact terms of Fe and V hyperfme fields showed 
a crucial role of valence contributions, being strongly de- 
pendent on local magnetic moments arrangements. The 
fact that only simple magnetic structures were accounted 
for the KKR calculations, did not allow to expect fully 
realistic determination of H va i . The enhancement of the 
H va i contribution which either lowers \Hhf\ on Fe atoms 
and increases this quantity on V atoms, would eventually 
lead to a better agreement with the experimental data. 
In spite of the aforementioned discrepancies between ex- 
perimental and calculated values, the almost linear de- 
pendences of < \Hp e \ > (a*), determined in cr-FeV phase 
for the FM and the APM models, remains in satisfying 
agreement with measured one. 
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